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Abstract 

An early history of metal enrichment in the Galactic halo is studied. We investigate chemical inhomo- 
geneity by using a stochastic chemical evolution model. The model confronts with metallicity distribution 
function of long-lived halo stars which is found to be a clue to obtain the best model prescriptions. We 
find that the star formation in the halo virtually terminated by ~ 1 Gyr and that the halo has never 
been chemically homogeneous in its star formation history. This conclusion does not depend whether mass 
loss from the halo is taken into account or not. Observed ratios of a-elements with respect to iron do 
not show scatters on a [a/FeJ-pe/H] plane, but this does not imply that interstellar matter in the halo 
was homogeneous because a chemical evolution path on this diagram is degenerate in the star formation 
rate. On the other hand, apparent spread of [Sr/Fe] ratio among metal-poor halo stars does not reflect an 
inhomogeneous metal enrichment, instead it is due to a sharp increase in a production rate of strontium 
that is probably synthesised in slightly less massive stars than progenitor of iron-producing SN II. 



, 1. Introduction 

At very beginning of enrichment history of the Galac- 
tic halo, interstellar medium (ISM) was certainly chem- 
, ically inhomogeneous, because star forming regions at 
that time were locally confined and stellar ejecta did not 
spread uniformly throughout the halo. Only a limited 

■ number of Type II supernovae (SNII) exploded and a 
, mean spatial separation of SNII was much larger than 

■ a typical radius of single supernova remnant (Audouze 
[ & Silk 1995). Thus, the ISM was enriched locally. A 
' trace of such inhomogeneous ISM should be frozen in 
, extremely metal-deficient stars which we observe today 

■ (e.g. McWilliam 1997). 

Nevertheless, no clear evidence which suggests that 
the halo was inhomogeneous in the past has so far been 
found. McWiUiam et al. (1995), Audouze & Silk (1995), 
and McWilliam (1997) claimed that large scatters of 
heavy metal abundances (e.g., Sr, Y, Zr, Ba, and Eu) 
of extremely metal-poor halo stars (Beers, Preston, & 
Shectman 1992; McWiUiam et al. 1995; Ryan, Norris, 
& Beers 1996) imply that the early halo was inhomoge- 
neous. However, these scatters may simply reflect vari- 
ous lifetime of stars producing heavy metals, whose ori- 
gins are poorly understood yet. The abundance of each 
heavy metal should sharply rise at a time when bulk of 
massive stars that produce it start to explode. Thus stars 
formed during such epoch should show a large scatter in 
relative abundance of heavy metals with respect to iron 
while they have nearly identical iron content. Therefore 
large scatters of heavy metals at certain values of [Fe/H] 
do not necessarily imply an inhomogeneous ISM in an 



early halo. 

On the other hand, even if relative abundance of a- 
elements (e.g., O, Mg, Si, Ca, etc) with respect to iron 
do show little scatter at any value of [Fe/H], it does not 
necessarily mean that an early halo was homogeneous. 
Even if a star formation history of each region in the 
halo is different from each other, evolutionary paths on 
the [Q;/Fe]-[Fe/H] plane should be nearly identical in an 
early evolutionary stage during which SNII were dom- 
inant sources of nuclear enrichment. The evolutionary 
path depends on a star formation rate (SFR) very weak. 

Although little is known about an early halo, there 
seems to be a clear evidence for the inhomogeneous ISM 
in a solar neighbourhood disc; long-lived dwarfs (e.g., 
Edvardsson et al 1993) show a spread of ~ 1 dex in iron 
abundance along poorly defined age-metallicity relation. 
By relaxing a usual assumption of well mixing, Copi 
(1997) estimated chemical evolution of the solar neigh- 
bourhood in a Monte Carlo fashion. The solar neighbour- 
hood is modeled by 1000 independent regions. The evo- 
lutionary history of one particular region is determined 
randomly based on the SFR and the initial mass function 
(IMF). Although he successfully reproduced the observed 
spreads of elemental ratios, he failed to explain the scat- 
ter appeared in the age-metallicity relation. Bateman & 
Larson (1993) suggested that random walk processes of 
atomic and molecular gas clouds are dominant mixing 
mechanisms of iron in the present day solar neighbour- 
hood disc. Wielen, Fuchs, & Dettbarn (1996) suggested 
that stellar orbital diffusion in combination with radial 
abundance gradients can induce inhomogeneity in the 
ISM. However, van den Hoek & de Jong (1997) showed 
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that the stellar orbital diffusion does not explain the 

abundance variation sufRcicntly. They suggested instead 
that a sequential star formation and an infall of metal- 
deficient gas play an important role in preventing the 
ISM from being mixed. Wilmes & Koppen (1995) studied 
chemical evolution of isolated individual ISM parcels and 
showed that mixing is inefficient in the galactic disc. Pi- 
lyugin (1996) suggested that major galaxy mergers form 
multiple stellar populations of different metallicities in 
the disc. While these mechanisms can well explain the 
scatters along the age-metallicity relation, one should 
also keep in mind that stellar abundances may not re- 
flect abundances of ISM from which stars formed, be- 
cause chemical condensation processes such as grain for- 
mation in circumstellar envelope (e.g. Henning & Giirler 
1986) and/or thermal diffusion in stellar atmosphere (e.g. 
Bahcall & Pinsonncault 1996) might work efficiently. 

If the early Galactic halo was chemically inhomoge- 
neous, the metallicity spread should appear among halo 
stars of the same age. Unfortunately, the present day 
stellar isochrone fitting cannot distinguish age of stars 
born in a young halo. Therefore, it is not clear at all if 
the metallicity spread exists among coeval halo stars. In- 
stead, we show in this paper that one can use an observed 
cumulative metallicity distribution function of halo stars 
to constrain stochastic chemical evolution models for an 
early halo and indicate that the Galactic halo has never 
been chemically homogeneous in its history. 

In section 2, we give prescriptions for our stochastic 
chemical evolution models, and in section 3 we describe 
behaviours of theoretical metallicity distributions. We 
confront model results with observational data in section 
4 and give discussions and conclusions in section 5 and 
6, respectively. 

2. Model 

2.1. Outline 

We have built up stochastic chemical evolution mod- 
els for the Galactic halo. Wc assume that the halo is 
spheroidal and that the gas of primordial abundance dis- 
tributed uniformly at the beginning. We divide; the; halo 
into many cubic blocks having a volume of = 47ri?;^/3, 
where Rm is a radius of supershell caused by SN II that 
contributed dominantly to an early stage of chemical evo- 
lution. Progenitor of SN II are massive stars which tend 
to form in star clusters and associations (Blaauw 1964; 
Humphreys 1978; Heiles 1987). Supernova remnants in 
OB associations and star clusters are suggested to pro- 
duce surrounding shell structures, or supershells (e.g., 
Tenorio-Tangle & Bodenheimer 1988). Cash et al. (1980) 
predicted that the supershell is produced after a series of 
SN II explosions. A radius of supershell can roughly be 
given as a maximum size of the area enriched by syn- 
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thesised heavy elements. Thus we assume that a volume 

of the block is equal to that of the supershell. The ISM 
in the supershell should be well-mixed, because stellar 
winds and SNII explosions induce the Rayleigh- Taylor 
and the Kelvin-Helmholtz instabilities (Tenorio-Taglc & 
Bodenheimer 1988; Allen & Burton 1993). The time scale 
of mixing is less than 0.0015 Gyr in a case of isonised gas 
(Roy & Kunth 1995), thus is much shorter than a typi- 
cal lifetime (~ 0.01 Gyr) of the OB association. Once the 
mixing took place, the ISM in a block uniformly enriched 
by newly synthesised elements. Therefore, individual re- 
gions can be regarded as independent one-zones and a 
standard chemical evolution model can be applied for 
each of them. 

Before an onset of star formation, all blocks have iden- 
tical physical properties. We therefore follow chemical 
enrichment histories of 1000 neighbouring blocks instead 
of the all. Unless a number of blocks considered is too 
small, the resulting behaviour of each block is sufficiently 
stable. By simulating stochastic chemical evolution for 
1000 blocks simultaneously, we calculate time variations 
of abundances in individual blocks and study chemical 
inhomogeneity in the halo. 

2.2. Block Size 

In an analogy to an expansion of single supernova, we 

roughly estimate a radius of the supershell as follows, 
although more detailed models for the supershells were 
published by several authors (Bruhweiler et al. 1980; Mc- 
Cray & Kafatos 1987; Mac Low & McCray 1988). Cioffi, 
McKee, & Bertschinger (1988) gave an analytical esti- 
mate for a radius Vm of a supernova remnant (SNR): 

= 69 E°,'^ v^r' (^/^o)"" "' PC, (1) 

where E^i is an initial kinetic energy of the SNR in units 
of 10^^ erg, n is the number density of surrounding ISM, 
wio is an ISM velocity dispersion in units of 10 km s"'^, 
and Z/Zq is the ISM heavy element abundance. We cal- 
culate a radius Rm of the supershell by replacing E^i 
with nil Ec,i^ where nu is the number of SNII in an OB 
association. Hereafter we adopt E^i = 1 (Shigeyama, 
Nomoto, & Hashimoto 1988). We assume that the ve- 
locity dispersion of ISM is equal to the sound velocity, 
and adopt ~ 10 kms~^ or vio — 1, since the tempera- 
ture should be near 10^ K (Hoyle 1953; Silk 1977) in the 
initial halo. Since we are interested in chemical evolu- 
tion in the halo, we fix Z/Zq = 0.06 which corresponds 
to [0/Fe]~ 0.4 (Barbuy 1988; Nissen et al. 1994) and 
[Fe/H]= —1.6; i.e.. a peak iron abundance of the metal- 
licity distribution function obtained for long-lived stars 
in the halo (Laird et al. 1988). 

Stars having M > 8Mq, corresponding to main- 
sequence spectral type B3, eventually become SNII. The 
number of supernovae nu is the number of such stars 
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per OB association. We adopt nu — 40 according to 
Heiles (1987) who derived the average number of stars 
with M > 8Mq per clusters from actual counting of O 
stars in clusters in the solar neighbourhood together with 
the IMF derived from clusters in the solar vicinity. The 
number Heiles (1987) got was ^ 28, which was then cor- 
rected for runaway O stars that will deposit energy within 
the supershell. Significant deviations from the average 
value of nil are apparently rare although nu may cover 
a wide range of values (Tenorio- Tangle & Bodenheimer 
1988). 

The number density of the ISM is given as n 
Mil/ (fimuVh), where mn is the proton mass, Mh is the 
mass, Vh — ^T^Rh halo volume, and /i = 1.3 is a 

mean molecular weight corresponding to the primordial 
compositions. 

Saito (1979) derived an empirical relation between the 
binding energy ^Iq and the mass Mq by analysing surface 
brightness distributions and line-of-sight velocity disper- 
sions of spheroidal systems; 



= 1.66 X 10^ 



Mg 



-I 1.45 



IQi^Mp 



erg. 



(2) 



Under the assumption of spherical geometry, the radius 
R{Mq) of a galaxy is given as 



R{Mg) = 26.1 



Mr 



1012 Mr 



0.55 



kpc. 



(3) 



We assume that the early Galactic halo follows the same 
mass-radius relation given by equation (H). The virial 
theorem tells that a system with initially no kinetic en- 
ergy attains virial equilibrium by reducing a radius to 
half the initial value. Since we study chemical evolution 
at the very beginning, the radius of halo i?h should be 
taken as 2i?(Mh). 

If we assume the halo was Mh = 4 • 10^^ M© (Fich 
& Tremaine 1991), then we obtain i?h = 2i?(Mh) = 
31.5 kpc, and n = O.lcm^'^. According to equation (1), 
the supershell radius i?ni depends weakly on n. Since 
the amount of gas condensed into stars in the early stage 
of halo evolution is at most 25 percent, i?m changes 
only by 12 percent. We therefore assume n constant in 
time. Putting nuE^i = 40, n = O.lcm"'^, viq ~ 1, and 
Z/Zq = 0.06 into equation (1), we obtain i?m ^ 660 pc 
and l\y = {AnR^/Sy^^ = 1 kpc for our standard model. 

A time tm when a SNR merges with the surrounding 
ISM is given by (Ciofh et al. 1988), 



;ro^-« (Z/Zo)-°-05 yr, (4) 



and a cooling time tc inside a SNR (Cox 1972) as, 
U - 5.7 • 10* E*^i" 71-9/1^ yr. 



(5) 



A lifetime of supershell and a cooling time Tc in- 
side the supershell can be derived by replacing E^i in 



equations (4) and (5) with nuE^i, respectively. As- 
suming the same parameters discussed above, we obtain 
~ 1.7 • 10-2 Gyr and ~ 4.6 • 10"^ Gyr. We note 
that is much longer than T^. 

2.3. Star Formation Probability 

Considering an idealised situation that the OB associ- 
ations distribute uniformly in the halo, we assume that 
a star formation probability Psf {t) in a block for a time 
interval t and t + is given as. 



PsY{t) 



1J_ 

Vh nil 



t+At 



Nii{t')dt' 



An \RhJ nil Jt 



Nii{t')dt', (6) 



where Nii{t) is the total number of progenitor of SNII in 
the halo and is given as. 



Nu{t) =. Mi,C{t)hi, 



(7) 



where C{t) is the SFR per unit mass and /// is the num- 
ber fraction of SNII defined as. 



/ii 



{m)m ^dm, 



(8) 



with mi =0.1 Mq, niu — 50 Mq, and mn — 10 Mq as 
the lower and the upper mass limits, and the lower SNII 
mass limit, respectively. The IMF adopted here has mass 
spectrum independent of time: 

0(™) = ^\^\r-i (mi <m<mu). (9) 

We adopt the Salpeter IMF (Salpeter 1955) which has a 
slope X = 1.35 in this definition. For C{t), we adopt the 
Schmidt law (Schmidt 1959): 



c{t) = G{ty {p = i), 



(10) 



where to and G{t) are time scale of star formation and 
a gas fraction in the halo, respectively. Psf(0 thus 
finally given as: 



(11) 



We adopt mu — lOAf© in equation (8), which is 
slightly larger than mu — 8Mq which Heiles (1987) 
adopted. Thus, equation (11) may underestimate PgF(t). 
A difference corresponding to the adopted mu is ~ 12 
percent. However, an exact value of mu is still uncertain 
and, therefore, we adopt mu = lOM© for our standard 
model and change it later as a free parameter. 
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2.4- Chemical Evolution 

For the blocks in which stars are born (hereafter the 
star forming blocks), wc apply a chemical evolution 
model developed by Pagel & Tautvaisiene (1995: here- 
after PT95). We adopt the iron as a tracer of chemical 
evolution. In PT95, the iron abundance is calculated 
by using an instantaneous recycling approximation and 
a delayed production approximation which assumes that 
the iron is additionally produced by Type la supernovae 
(SNIa) after a fixed time lag ri. We assume that the 
SFR c{t) in a block is proportional to the gas fraction 
g{t). The iron abundance .^Fe(i) is given by, 

dt 

^( -Zp^{t)c{t)+poc{t}, (t<Tl), 
\ -Zpe(t)c(t)+poc(t)+Plc(t-Tl), (t>Tl), 

where po and pi are yields of iron corresponding to 
the instantaneous recycling approximation and the de- 
layed production approximation, respectively. The first 
term on the right hand side (r.h.s.) represents the net 
amount of iron locked into newly formed stars. The sec- 
ond and the third terms are corresponding to the instan- 
taneous recycling approximation and the delayed produc- 
tion approximation, respectively. We adopt po = 0.28, 
Pi = 0.42, and ti = 1.3 Gyr according to PT95. Since 
the theoretical yields of the iron are uncertain due to dif- 
ficulties in modeling the explosion mechanism (Timmes, 
Woosley, & Weaver 1995), we use the yields of PT95 
which are calibrated empirically. 

The SFR in each block c{t) is determined by an as- 
sumed number of SN II exploded in a block during a time 
interval t and t + At: 

/t+At 
rubfiicit') dt' = nil, (13) 

where mb = Mh/^/Vh — 3 • 10^ Mq is the initial gas mass 
of each block. 

2.5. Standard Models 

Table 1 gives a list of parameters adopted for our stan- 
dard models SI, S2, and S3. In model SI we assume 
that the OB associations are formed randomly at every 
At ^ = 0.006 Gyr. In model 82 we introduce an 
effect of periodic mixing due to the turbulent diffusion. 
A particle with a sound velocity v = 10km s~^ would 
across a block of size Zb = 1 kpc in = 0.1 Gyr. Thus, 
we assume in this model that every 0.1 Gyr the turbu- 
lent diffusion mixes the 27 blocks surrounding the star 
forming block periodically. 

The supershell structures are often mentioned as a trig- 
ger of massive star formation. The evidences have been 
discussed by many authors (Blaauw 1964; Elmegreen & 
Lada 1977; Lada, Blitz, and Elmegreen 1979; Elmegreen 
1982; 1985a, b). If OB associations form in this way, 
they should be born only in the surrounding area of the 



Table 1. Parameters of standard model. 





4 10" Mq 


mass of the halo 


Rh 


31.5 kpc 


radius of the halo 


LO 


5 Gyr 


time scale of star formation in the halo 


lb 


1 kpc 


size of a block 


VlO 


1 


velocity dispersion of ISM 






in units of 10 km 


mi 


O.lAfc? 


IMF lower mass limit 


mil 


IOMq 


lower mass limit for SNII 


nil 


40 


number of SNe II in a block 


At 


6 • 10® yr 


duration of star formation 



star forming regions and chemical enrichment should be 
locally confined there. The effects of this stimulated star 
formation are studied in model S3. We first distribute OB 
associations randomly and then assign randomly a star 
forming block for the next generation among 27 blocks 
adjoining to each star forming block of the first genera- 
tion and so on. In this way, the star formation propagates 
from one block to another. 

The time step adopted for each simulation is 0.001 Gyr. 
Calculations are stopped at 1 Gyr, since our main interest 
is in early chemical evolution of the halo. 

The time scale of star formation cj in the halo is quite 
uncertain. On the contrary, it is quite often assumed that 
(jj in elliptical galaxies is of the order of free-fall time 
ts = (37r/32Gp)i/2 (Larson 1974). One could assume 
that u! in the halo is the same as that in giant ellipticals. 
However, mean stellar metallicities of giant ellipticals are 
nearly solar (Arimoto et al. 1997), while that of Galactic 
halo stars is ^ 1/40 solar (Ryan & Norris 1991; Carney 
et al. 1996). Both in giant ellipticals and in the galactic 
halo, Mg seems to be enhanced with respect to Fe by 
at least a factor of 2 (McWilliam et al 1995; Ryan et al 
1996). This implies that star formation in both systems 
is characterised by the same IMF and that the star for- 
mation stopped before the onset of SNIa explosions (~ 1 
Gyr). If this is the case, lu in the Galactic halo should 
be 40 times longer than that of giant ellipticals. Thus 
we obtain uj ~ 40tff ^ 6 Gyr for our halo model with 
Mh = 4 • 10" Mq and i?h = 31.5 kpc. Since lo - 6 Gyr 
is vciry close to ^ 5 Gyr derived for the solar neigh- 
bourhood disc (Arimoto, Yoshii, & Takahara 1992), we 
shall adopt w = 5 Gyr in our standard models. In other 
words, we adopt the same SFR fo the Galactic halo and 
the solar neighbourhood disc. 

3. Metallicity Distribution 

Figure 1 shows the number fraction of blocks which 
at least once have experienced the chemical enrichment. 
Model 81 predicts that all the blocks have been enriched 
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Fig. 1.. The evolution of the number fraction of blocks which at 
least once have experienced chemical enrichment in the stan- 
dard models SI, S2, and S3. 

before ~ 0.2 Gyr. In model S2, this happens earlier than 
model SI, because the iron spreads beyond the super- 
shells due to the turbulent diffusion, thus much wider 
area are enriched even if the SFRs are the same as in 
model SI. On the other hand, nearly 8% of the block in 
model S3 has never been enriched till 1 Gyr, this is be- 
cause the stimulated star formation tends to localise the 
chemical enrichment. 

Figures 2a and 2b show the evolution of mean and me- 
dian iron abundances, respectively. Hereafter we discuss 
statistical properties except for the first 0.04 Gyr, during 
which the number of enriched blocks is too small to define 
the statistical properties. Figures 2a and 2b indicate that 
the median metalHcities of models SI and S2 arc nearly 
the same as the mean metallicities. This suggests that 
the metallicity distributions of the two models have sym- 
metric shapes. The median of model S3 is always higher 
than the mean, indicating that the metallicity distribu- 
tion has a tail toward lower metallicity. 

Figure 3 illustrates a frequency distribution of iron 
abundance of the 1000 blocks, iV([Fo/H]), of model SI 
at 0.01, 0.05, 0.1, 0.2, 0.5, and 1 Gyr. For an illustrating 
purpose, N{—5 <[Fe/H]< —4.8) is artificially assigned 
for a fraction of metal-free blocks. Figure 3 shows that 
the metallicity distribution of model SI has roughly 1.3 
dex spread of the iron abundance, suggesting that the 
ISM in the halo was inhomogeneous till 1 Gyr. 

Figure 4 gives the evolution of standard deviation of 
the metallicity distribution obtained by the standard 
models. In model SI, the standard deviation keeps 
nearly constant at s{t) ~ 0.24 after ~ 0.2 Gyr, which 
corresponds to roughly 1.3 dex dispersion of the iron 
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(a) 


The evolution of the mean iron abundance in the 



standard models SI, S2, and S3, (b) The same as figure 2(a), 
but for the median iron abundance. 



abundance. Model S3 suggests much stronger abun- 
dance spread than model SI. The standard deviation of 
model S3 increases with time and exceeds 0.5 at 1 Gyr, 
which is roughly equivalent to 2 dex dispersion of the 
iron abundance among all the enriched blocks. On the 
other hand, the standard deviation of model S2 decreases 
quickly to nearly zero, showing that the ISM was homoge- 
neous from the very beginning of its chemical evolution. 



Now we focus on the shape of the metallicity distribu- 
tion. Figures 5 a and 5 b illustrate the skewness and the 
kurtosis of the metallicity distributions of the standard 
models, respectively (e.g. Stuart & Ord 1987). The skew- 
ness characterises a degree of asymmetry of a distribution 
around its mean. It characterises only a shape of the dis- 
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Fig. 3.. The frequency distribution of iron abundance of the 1000 blocks in model SI at 0.01, 0.05, 0.1, 0.2, 0.5, and 1 Gyr. 



tribution. The definition is 

1 ^ - " 

Skew{xi, X2,...., xn) = j^^i^^-j^f, (14) 

i—l 

where x and s are the mean and the standard deviation 
of the measured value Xi, X2,---, x^. A positive value of 
skewness signifies a distribution with an asymmetric tail 
extending out towards larger x, while a negative value 
signifies a distribution whose tail extends out towards 
smaller x. The definition of the kurtosis is 

Kurt(xi, X2, xr,) = {- V(^^)n - 3, (15) 

iV f— f s 

where the first term of r.h.s. becomes 3 for a Gaussian 

distribution. The kurtosis measures the relative peakcd- 
ness or flatness of a distribution relative to the Gaussian 
distribution. The distribution with positive kurtosis has 



the outline of the Matterhorn for example. The distri- 
bution with negative kurtosis is outlined of a lump of 
meat-loaf. Figure 5a gives the evolution of skewness of 
the standard models. The skewness of model SI is nearly 
zero, which means a symmetric shape of the mctallicity 
distribution. The skewness of model S2 behaves irregu- 
larly due to nearly null standard deviation. The skew- 
ness of model S3 is always negative and decreases as times 
goes on, showing a tail of the metallicity distribution ex- 
tending toward lower metallicity at later stages. This is 
due to the localised chemical enrichment caused by the 
stimulated star formation. Figure 5b shows the evolution 
of kurtosis. The kurtosis of model SI is nearly equal to 
zero, thus the metallicity distribution has a Gaussian like 
shape. On the contrary, the kurtosis of model S2 is very 
large. This model suggests the homogeneous ISM. since 
the kurtosis should be infinitely large if samples have the 
same value. In model S3, the kurtosis gradually increases 
up to ~ 2, which means that the metallicity distribution 
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models SI and S3 show that the ISM in the halo was 

not well-mixed till at least 1 Gyr. Indeed, the scatters 
of the iron abundance are roughly 1.4 dex in model SI 
and 2 dex in model S3, respectively. We have repeated 
the simulations several times and confirm that the model 
results show always the same tendencies. 

4. Observational Constraints 

Long-lived stars in the halo should keep the original 
metallicities of the ISM from which they formed. It is 
true that little observational information is available for 
constraining chemical evolution models of the halo, but 
we will show that a cumulative metallicity distribution 
function of the long-lived halo stars can potentially give 
a powerful clue to obtain the best model prescriptions. 



Fig. 4. . The evolution of standcird deviation of metallicity distri- 
bution in the standard models SI, S2, and S3. 



5 - I 






i (Gyr) 

Fig. 5.. (a) The evolution of skewness of metallicity distribu- 
tion of the standard models SI, S2 and S3, (b) The same as 
figure 5a, but for the kurtosis of metallicity distribution. 



has a sharp peak. Thus the skewness and the kurtosis of 
model S3 indicate that the metallicity distribution has a 
peaked shape with a tail extending to the lower metallic- 
ity. 

As a summary, only model S2 shows that the ISM was 
chemically homogeneous in the early halo, while both 



4-1- Comparison with the Standard Models 

Figures 6 and 7 present generalised metallicity distri- 
bution functions (GMDFs), P{[Fc/U]), of the standard 
models SI - S3 together with the empirical GMDFs of 
long-lived halo stars taken from Ryan & Norris (1991) 
and Carney et al. (1996). The definition of GMDF is 
given by Laird et al. (1988): 



Pix) 



1 



N 



^exp 



(x - Xj)"^ 

2^ 



(16) 



where a reads as a typical error in the observed values 

Xi, X2, Xjv. 

Quite often, the empirical metallicity distribution func- 
tion is presented in the form of a histogram. However, 
Searle & Zinn (1978) suggested that the bins of a con- 
ventional histogram should be replaced by a continuous 
distribution function such as a Gaussian to obtain a bet- 
ter approximation to the actual abundance distribution 
function, since the binning distorts the data. Therefore, 
we convolve the observed histograms of stellar metallicity 
distribution to include the uncertainty of observational 
data properly (Laird et al 1988) and convolve the the- 
oretical GMDFs by using equation (16) with a = 0.15 
(Ryan & Norris 1991; Carney et al. 1996). We note that 
the GMDFs of Ryan & Norris (1991) and Carney et al. 
(1996) are almost identical, except that the latter shows 
a small hump at [Fe/H] ~ —2.7. 

Figure 6 shows that model SI gives good fits to the ob- 
served GMDFs at 0.5 Gyr and 1 Gyr. The star formation 
in the halo must have virtually stopped at around 0.5 
Gyr, otherwise the resulting GMDF gives too high peak 
metallicity. This would happen if the gas escapes from 
the halo and accretes on to the bulge and disc dissipation- 
ally. A precise value of the epoch when the star formation 
terminated is of course model dependent. However, the 
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Fig. 6.. The generalised metallicity distribution functions (GMDFs) of model SI at 0.01, 0.05, 0.1, 0.2, 0.5, and 1 Gyr (thick solid 
line). Thin solid and dotted lines illustrate the empirical halo GMDFs taJjen from Ryan & Norris (1991) and Carney et al. (1996), 
respectively. 



GMDFs of halo stars strongly suggest that the star for- 
mation in the halo did not last longer than, say, 0.5-1 
Gyr. 

Upper and lower panels in figure 7 illustrate the 
GMDFs of models S2 and S3, respectively. Clearly, 
model S2 is inconsistent with the observations. Model S2 
gives too sharp GMDFs to reproduce the observations. 
We have also studied several alternative cases in which 
wc assume much longer time interval for the periodic 
turbulent mixing, much smaller size of the smoothed 
out area, and lower or higher probability of star forma- 
tion, and have confirmed that these models give much 
sharper GMDFs than the empirical ones; thus inconsis- 
tent with the observations. On the contrary, the GMDF 
of model S3 is consistent with those observed. Equiva- 
lent models to model S3, but with higher or lower SFRs, 
are also calculated. These models predict similar inho- 



mogeneous ISM enrichment. Their GMDFs are also con- 
sistent with the empirical ones at 0.1 Gyr (higher SFR) 
and 0.5 Gyr (lower SFR). The standard deviations of 
metallicity distribution function of these models increase 
with time, showing the same trends as those of model S3. 
Without additional informations, it is rather difficult to 
conclude which star formation mechanism, the sponta- 
neous star formation (model SI) or the stimulated one 
(model S3), is responsible for the early enrichment of the 
Galactic halo. 

To judge a goodness of the model fit to the empirical 
GMDFs, we have performed the statistics. Let Xi 
and Ui be the theoretical and the observed P([Fe/H]i), 
respectively. Then, the value is defined as. 



E 



jxi - Vif 
Xi+Vi 



(17) 
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Fig. 7.. The GMDFs of model S2 (upper three panels) cind model S3 (lower three panels) at 0.1, 0.5, and 1 Gyr. Thick solid lines show 
the theoretical GMDFs. Thin solid and dotted lines have the same meaning as in Fig. 6. 



In columns (8)-(10) of tabic 2, wc give the values of the 
standard models at 0.1, 0.5, and 1 Gyr, respectively. The 
best fit is realised by model SI at 0.5 Gyr, but model SI 
at 1 Gyr and model S3 at 0.5 Gyr also give a reasonable 
fit. 

In addition to the standard models, wc show additional 
7 models A - G to demonstrate the parameter depen- 
dence of our results. All the 7 models are equivalent 
to model SI, but with different parameters, i.e., these 
models assume spontaneous star formation and no tur- 
bulent mixing. The parameters of these 7 models are 
summarised in tabic 2. In this table, column (1) gives an 
identification of the model, columns (2) and (3) give the 
adopted values of Ih/Rh and u}~^At, respectively. Col- 
umn (4) gives the parameter that is different from the 
standard models. Columns (5) - (7) indicate the stan- 
dard deviation s{t) at 0.1, 0.5, and IGyr, respectively. 
The values for the best fit epochs of the standard mod- 
els are 11.01 (model SI) and 11.89 (model S3) at 0.5 Gyr. 
Therefore, we consider a good fit is realised if the model 



gives the less than 12 in columns (8) - (10). 

In models A - C, we study how the evolution of s{t) and 

X^ depend on Ih/Rh and oj~^At. Two cases of Ih/Rh = 

0. 01 and 0.03 are studied. For LU~^At, we consider three 
cases, i.e. w-^At = 0.6- 10"^, 1.2-10"^, and 3-10-2. The 
standard deviation of model A is always larger than that 
of model SI, which suggests that the ISM is always chem- 
ically inhomogeneous. In model B, higher SFRs lessen 
the differences of chemical enrichment among the blocks. 
The standard deviation of model B is nearly equal to 0.1 
at 1 Gyr, predicting relatively well-mixed ISM with re- 
spect to model SI. However, the of model B indicates 
that this model is inconsistent with the observations. In 
model C, given smaller Ih/Rh, the best fit epoch to the 
observations comes earlier than model SI. The gas in 
the star forming blocks of model C is converted into stars 
more efficiently than model SI, since equation (14) indi- 
cates that the SFR in a block is proportional to nn/wb. 

1. e., {Ih/Rh)^"^- In modelD, we reduce the number of 
SNII progenitor in a single OB association and assume 
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Table 2. Models and results. 
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V- ; 
0..J 
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0.1 


().") 


I 


SI 


0.03 


1.2 




0.27 


0.24 


0.24 


75.11 


11.01 


11.85 


S2 


0.03 


1.2 




0.27 


0.01 


0.01 


76.82 


26.78 


22.63 


S3 


0.03 


1.2 




0.33 


0.45 


0.54 


68.16 


11.89 


20.36 


A 


0.03 


0.6 




0.34 


0.34 


0.29 


53.52 


4.68 


15.82 


B 


0.03 


3 




0.18 


0.10 


0.10 


95.12 


24.40 


16.79 


C 


0.01 


1.2 




0.32 


0.44 


0.42 


9.25 


11.68 


19.73 


D 


0.03 


1.2 


nil = 20 


0.15 


0.11 


0.11 


92.31 


21.39 


17.19 


E 


0.03 


1.2 


mii=8M0 


0.24 


0.16 


0.17 


83.33 


16.70 


14.02 


F 


0.03 


1.2 


mi = O.O5M0 


0.30 


0.33 


0.30 


62.80 


6.19 


11.77 


G 


0.03 


1.2 


Mh = 2 • 10" Mg 


0.30 


0.42 


0.37 


49.38 


4.00 


13.06 



Model S2 assumes propagation of star formation. Model S3 assumes the mixing of ISM by turbulent diffusion whose priod is 
taken to be 10* yr. See text in detail. 



nil = 20 instead of 40 (Blaauw 1964). This increases 

the probability of star formation (sec equation (11)) and 
homogenises the ISM quickly. The of model D in- 
dicates poor fits to the observations. Model E assumes 
mil = 8Mq. The number fraction of SNII progenitor 
in model E is increased and the SFR becomes twice of 
model SI. Thus, both the evolutionary behaviour of stan- 
dard deviation and the fits arc very similar to those 
of model B. We can reject this model, since the gives 
very poor fits. Models F and G assume mi = O.O5M0 and 
Mh = 2 • 10^"'^ Mq, respectively. Equation (11) indicates 
that the probability of star formation P{t) in a block, 
namely the SFR in the halo, is proportional to the num- 
ber fraction /n of SN II progenitor and the mass of the 
halo. Thus in models F and G with the lower SFR in the 
halo, the standard deviation of metallicity distribution 
function becomes slightly larger than that of model SI. 
The x^ values of models F and G show that these models 
give a good fit to the observed GMDFs in the halo at 
0.5 Gyr. 

In summary, the different parameters change the SFR 
in the halo. If the SFR is low, the standard deviation 
of metallicity distribution becomes large, i.e. the ISM is 
inhomogeneous, and good fits to the observed GMDFs 
are resulted. While if the SFR in the halo is high, the 
standard deviation becomes small (well-mixed ISM) and 
the resulting GMDFs show poor fits to the observations, 
since stars born in the similarly enriched ISM dominate 
the cumulative metallicity distribution function at later 
phase. Thus the predicted GMDF becomes too narrow 
to be consistent with the observations. 

4-2. Comparison with Mass-Loss Models 

In previous sections, we assume that the halo was a 
closed system. Hartwick (1976) showed that a simple 



model predicts too many metal-poor stars and is incon- 
sistent with the cumulative distribution of the iron abun- 
dances for metal-poor ([Fe/H]< —0.52) globular clusters 
in the Galactic halo. He modified the simple model and 
proposed a model which allows for gas outflow from the 
halo. We therefore consider effects of mass loss from the 
halo in this subsection. The gas will eventually either es- 
cape to intergalacitc space or to accrete onto the galactic 
plane. 

Since the gas removal from the halo could be related to 

energy injection from hot OB stars and SNII (Hartwick 
1976), we assume that the mass loss rate dD/dt is pro- 
portional to the SFR C{t): 



dP 
~dt 



bC{t). 



(18) 



The mass loss rate h and the SFR are taken from 

Hartwick (1976), i.e., 6 = 10 and = 0.3, respectively. 
We assume the same gas outflow rate for all the blocks in 
the halo. In model MLl, the parameters are the same as 
those of model SI except for b and uj~'^. We assume that 
OB associations are born randomly in space and do not 
consider the ISM mixing by turbulent diffusion. While 
in model ML2, we study the effects of turbulent mixing. 
We assume periodic ISM mixing, similarly to model S2. 
The period is taken as 0.1 Gyr, the same as that of model 
S2. Other prescriptions for model ML2 are same as those 
of model MLl. 

Figures 8a and 8b show the evolution of mean iron 
abundance and that of standard deviation of the mass 
loss models, respectively. Figure 9 shows the theoreti- 
cal GMDFs plotted with the observed GMDFs. These 
figures show that the GMDFs obtained by model MLl 
agree with the observed GMDFs and the ISM is always 
inhomogeneous in this model. On the other hand, the 



10 C. Ikuta & N. Arimoto 




i (Gyr) 



Fig. 8.. (a) The same as figure 2(a), but for mass loss models 
MLl and ML2. The SFR u)^'^ = 0.3 and the mass-loss rate 
6 = 10 arc the same as those of Hartwick (1976). (b) The 
same as figure 4, but for mass-loss models MLl and ML2. 



GMDFs of model ML2 are too sharp and are inconsistent 
with the observations, because the ISM in this model is 
wcll-rnixcd. In model MLl, the best fit is achieved at an 
epoch later than model SI. In the case of one-zone model, 
a shift of the peak abundance can be explained by a de- 
crease of effective yield caused by the mass-loss (Hartwick 
1976; Pagel 1992). Although our model is not one-zone, 
the mass-loss gives a similar effect to the peak abundance 
of the cumulative metallicity distribution function. 

The GMDFs obtained by assuming & = 10 and lo~^ = 
0.2 (the same SFR as that of model SI) also fit to the 
observed GMDFs and the best fit is realised at around 
1 Gyr. The chemical enrichment in this model is also 
inhomogeneous. A model adopting 6=5 also gives con- 
sistent GMDFs with the observations, although the best 
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fit epoch appears earlier than model MLl due to a larger 
effective yield. Gas outflow may occtir in the vicinity of 
star forming regions, if major energy sources are massive 
stars. When we take this into account and consider a 
model in which gas is expelled only from the star form- 
ing blocks, we obtain the same results on a fitting good- 
ness and inhomogeneity of the ISM as those obtained by 
models MLl and ML2. 

As a conclusion, whether or not mass loss is taken into 
account, the ISM in the halo at the beginning must be 
inhomogeneous . 

5. Discussions 

Observations show neither clear age- metallicity rela- 
tion nor clear evidence for metallicity gradient in the 
halo. Because of these observational results. Galactic 
halo formation scenario favoured recently (e.g Carney et 
al. 1996) is the one proposed by Searle & Zinn (1978), i.e. 
accretion of small proto-galactic fragments contributed to 
the halo population. However, in the poorly mixed halo, 
the age-metallicity relation and the metallicity gradient 
cannot be expected even if the Galaxy formed from one 
massive cloud. Thus it is not necessary to introduce the 
accretion of proto-galactic fragments to interpret the lack 
of the agc-mctallicity relation and the metallicity gradi- 
ent in the halo. 

The scatters observed in the relative abundances 
of neutron capture elements to the iron, particularly 
[Sr/Fe], are often claimed as evidences for the inhomo- 
geneous enrichment of the ISM in the halo at the very 
beginning. In figure 10, Sr abundances of halo stars 
taken from Gratton & Sneden (1988; 1994), Magain 
(1989), McWilliam et al. (1995), and Ryan et al. (1996) 
are plotted. The observations show that a majority of 
the metal-poor stars with [Fe/H] < —2.5 distribute at 
— 1.5 < [Sr/Fe] < 0.6. The relative abundance of Sr tends 
to decrease clearly with an increase of the iron abun- 
dance. 

We show two theoretical evolutionary paths of Sr on 

figure 10. The models are calculated by adopting the 
same parameters as those of models SI (left path) and C 
(right path). The SFR of model SI is nearly one tenth 
of model C. Co-existence of regions with different SFR is 
a view of the chemical evolution in the halo considered 
here. We have simplified the situation and assumed the 
same SFR in each block. However, different SFR in each 
star forming region is more realistic. These models are 
consistent with the observed behaviour of [Sr/Fe] and the 
empirical metallicity distribution ftmctions. 

Following PT95 and Pagel & Tautvaisiene (1997; here- 
after PT97)), we calculate the chemical evolution of Sr: 

-rAgZs.) = -Zsr{t)c{t) +poc(i) +^ftc(t-Ti),(19) 
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Fig. 9.. The GMDFs of model MLl (upper three panels) and model ML2 (lower three panels) at 0.1, 0.5, and 1 Gyr. Thick solid lines 
show the theoretical GMDFs. Thin solid and dotted lines have the same meaning as in Fig. 6. 



where g{t), c{t), and ^Sr(i) are the gas fraction, the SFR, 
and abundance of Sr in a block, respectively. The first 
and the second terms on the r.h.s. have the same meaning 
as those in equation (12). The third term comes from the 
delayed production, where pi is the yield corresponding 
to the fixed time lag t^. Here po = 0.01, pi — 0.08, 
P2 = 0.39, P3 = 0.23, Ti = 0.023 Gyr, T2 = 0.025 Gyr, 
and T3 = 2.7 Gyr are assumed (PT97). The evolution of 
the iron abundance is calculated by equation (12). 

The theoretical paths in figure 10 show that [Sr/Fe] be- 
gins to increase rapidly at [Fe/H] ^ —3.4 in model SI and 
[Fe/H] ~ —2.3 in model C, depending on the SFR. These 
paths roughly outline the dispersion of observed [Sr/Fe] 
at [Fe/H] < —2.5. The models are calculated with a time 
step of 0.001 Gyr, which is shown by filled pentagons on 
the theoretical paths in figure 10. The chemical evolu- 
tion from [Sr/Fe] ~ —1.3 to [Sr/Fe] ~ should take place 
within a short time interval (~ 0.012 Gyr). After the 
rapid increase, the paths show temporal increase and de- 
crease, and then converge. The declines of theoretical 



[Sr/Fe] are not because of delayed production of the iron 
(in other words, contribution from SNIa which produces 
the bulk of iron in the solar neighbourhood), since we 
assume that the time delay of iron r in equation (|l^) is 
equal to 1.3 Gyr (see section 2.3) and consider chemical 
evolution before 1 Gyr. The declines of [Sr/Fe] are due 
to onsets of the next formation of OB associations. Thus 
different SFR in each star forming region can explain the 
trends defined by a majority of observed stars on [Sr/Fe] 
vs [Fe/H] diagram. 

Apparently the models are inconsistent with the ob- 
servational data on figure 10. McWillam et al (1995), 
Sneden et al (1996), and McWillam (1998) have shown 
that the heavy element abundances in the super-solar 
[Sr/Fe] stars are dominated by the r-process abundance 
patterns. The r-process elements can only be produced 
by Type II supernova events. Therefore one may argue 
that the time-delay model for Sr is inadequate to ex- 
plain the observations. If the strontium is produced in 
all progenitors of SNII, however, [Sr/Fe] values observed 
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209 ([Fe/H] 3.2), Primas, Molaro, & Castelli (1994) 

and McWilliam et al (1995) reported [Sr/Fe]= 0.92 and 
[Sr/Fe]= —0.07, respectively, although they used the 
same lines and had the same S/N. 



*Gratton & Sneden (1988) - 
OGratton & Sneden (1994) - 
OMagain (1989) 
□ McWilliam el ai. (1995) — 
A Ryan et al. ( 



- 




Fig. 10.. Chemical evolution of strontium. Observational data 
are taken from Gratton & Sneden (1988; 1994), Magaiii (1989), 
McWilliam ct al. (1995), and Ryan et al. (1996). Two theo- 
retical paths correspond to different SFR in a block. The left 
and right paths are obtained by adopting the same parameters 
as those of models SI and C, respectively. Each filled pentagon 
on the paths indicates a time step 0.001 Gyr. 



in metal-poor stars should be always super-solar. On 
the contrary, as we mentioned before, figure 10 shows 
that [Sr/Fe] tends to decrease at lower iron abundance 
([Fe/H] < —2.5). This drop suggests that the strontium 
was formed slightly later than the iron (see Mathews, 
Bazan, & Cowan 1992). 

If the large scatter of [Sr/Fe] observed for stars of iron 
abundance in the range of —3.6 <[Fe/H]< —2.6 is in- 
deed due to a sharp increase of the Sr production rate, 
the Sr abundances of four extremely metal-dificient stars 
with [Fe/H]~ —4 seem to show some evidence against 
it. However, the iron abundance measurement of such 
very low metal stars is extremely difficult and these iron 
abundances could well be as large as [Fe/H]> —3.6 (M. 
Spite, private communications). 

The theoretical [Sr/Fe] is too low to fit to the observed 
extremely high [Sr/Fe] value. However, the abundance of 
these stars might not reflect the composition of ISM from 
which they formed (e.g. McClure 1984). We should re- 
ject CH stars when we discuss the Galactic chemical evo- 
lution, because atomosphers of CH stars are thought to 
have been enriched by elements transfered from evolved 
companion AGB stars (McClure 1984). McWillam et al 
(1995) and McWilliam (1998) reported that CS22898- 
027 ([Fe/H]= -2.36) and CS22947-187 ([Fe/H]= -2.5) 
are probably CH stars. Observational errors might also 
reflect the [Sr/Fe]- [Fe/H] diagram. For a star CS22891- 



6. Conclusion 

A stochastic chemical evolution model has been built 

to study an early history of metal enrichment in the 
Galactic halo. The metallicity distribution function of 
long-lived halo stars is found to be a clue to obtain the 
best model prescriptions. Wc find that the star forma- 
tion in the halo virtually terminated by ~ 1 Gyr and that 
the halo has never been chemically homogeneous in its 
star formation history. The star formation in the halo 
could either be spontaneous or stimulated, which keep 
the halo always inhomogeneous and the turbulent mix- 
ing is found to be inefficient. This conclusion docs not 
depend whether the mass loss from the halo is taken into 
account or not. The observed ratios of the a-elements 
with respect to the iron do not show scatters on the 
[Q;/Fe]-[Fe/H] plane, but this does not imply that the 
ISM in the halo was homogeneous because the chemical 
evolution path on this diagram is degenerate in the SFR. 
On the other hand, the apparent spread of [Sr/Fe] ratio 
among metal-poor halo stars does not reflect an inhomo- 
geneous metal enrichment, instead it is due to a sharp 
increase in the production rate of strontium that is prob- 
ably synthesised in slightly less massive stars than the 
progenitor of iron-producing SN II. 
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